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Abstract 

Near to filling fraction v = 1, the quantum Hall ferromagnet contains 
multiple Skyrmion spin excitations. This multi-Skyrmion system has a 
tremendously rich quantum-critical structure. This is simplified when 
Skyrmions are pinned by disorder. We calculate the nuclear relaxation 
rate in this case and compare the result with experiment. We discus how 
such measurements may be used to further probe the quantum-critical 
structure of the multi-Skyrmion system. 



At exact filling of a single Landau level the quantized Hall state forms 
an almost perfect ferromagnet. This quantum Hall ferromagnet (QHF) has 
some novel features due to the phenomenology of the underlying quantized Hall 
state. Magnetic vortices, or Skyrmions, in the QHF carry quantized electrical 
charge [Q. These Skyrmions are stabilised by a chemical potential so that the 
ground state slightly away from filling fraction v = 1 contains a finite density of 
them. This original proposal of Sondhi et al. has been confirmed in a number 
of experiments (l], gj. 

The T — phase diagram of the multi-Skyrmion system has been thor- 
oughly investigated. In the absence of disorder, crystalline arrangements are 
expected H Q. For filling fractions very close to v = 1 a triangular lattice is 
formed with a transition to a square lattice as the deviation from v = 1 is in- 
creased. The statistical mechanics of the possible melting transitions has been 
considered^. At the highest Skyrmion densities, zero point fluctuations are 
expected to give rise to a quantum- melted state gj. 

Despite this wealth of study, a complete account of the experimental observ- 
ables has not been achieved. For example, nuclear magnetic resonance provides 
one of the clearest probes of the spin polarisation in the QHF [|| . Although it is 
understood in general terms how low-energy spin-fluctuations of the Skyrmion 
system may enhance the relaxation of nuclear spins, attempts to calculate re- 
laxation rates have been flawed^. The fundamental physics, missed in other 
considerations, is the quantum-critical nature of fluctuations of the Skyrmion 
lattice. One immediate consequence of this quantum-criticallity is that the 



limits of temperature/frequency— > and frequency/temperature— > are very 
different. Typically, experimental probes are at frequencies much less than tem- 
perature and the latter limit is appropriate. This means that zero-temperature 
calculations cannot model experiments correctly]^]. 

Let us consider these points a little further. In our analysis below, we will 
find an underlying gappless XY-model governing orientational fluctuations of 
the multi-Skyrmion system. If a spinwave expansion is attempted for a gapp- 
less magnet at or below its critical dimension, the occupation of low-frequency 
modes is found to diverge. The constraint, fixing the magnitude of the lo- 
cal spin, restricts this divergence. Interplay between divergence and constraint 
gives rise to a finite (temperature dependent) correlation length, £(T), beyond 
which correlations of the magnet decay exponentially. The dynamics of the 
critical magnet are very different on length scales greater or less than £. On 
lcngthscales less than £, the groundstatc is ordered (albeit in a quantum super- 
position of all possible orientations due to long wavelength spin fluctuations). 
Fluctuations with wavelength less than £(T) may, therefore, be described by 
a modified spinwave expansion. On lengthscales greater than £, the ground- 
state is disordered and fluctuations are overdamped. This quantum relaxational 
dynamics is a striking feature of quantum-critical systems and leads to interest- 
ing universalities Regimes of renormalized classical and quantum activated 
behaviour at low-temperature, cross over to universal behaviour in the high- 
temperature quantum-critical regime. 

The Skyrmion spin-configuration consists of a vortex-like arrangement of in- 
plane components of spin with the z-component reversed in the centre of the 
Skyrmion and gradually increasing to match the ferromagnetic background at 
infinity. At large distances, the spin distribution decays exponentially to the 
ferromagnetic background on a length scale determined by the ratio of spin 



stiffness to Zeeman energy) 10 . An individual Skyrmion may be characterized 
completely by its position (i.e. the point at which the spin points in the opposite 
direction to the ferromagnetic background) its size (i.e. the number of flipped 
spins) and the orientation of the in-plane components of spin. The equilibrium 
size of the Skyrmion is determined by a balance between its coulomb and Zeeman 
energies Q (In the presence of a disorder potential, the potential energy of the 
Skyrmion also enters this balancing act fll|| ) . 

Consider a ferromagnet with a dilute distribution of Skyrmions. The normal 
modes of this system are relatively easy to identify. Firstly, ferromagnetic spin- 
waves propagate in-between the Skyrmions. The spectrum of these is gapped 
by the Zeeman energy and will be ignored from now on. Positional fluctuations, 
or phonon modes, of the Skyrmions are gapless in a pure system, but gapped 
when the lattice is pinned by disorder. Finally, fluctuations in the in-plane 
orientation and size must be considered. These two types of fluctuation are in- 
timately connected; rotating a Skyrmion changes its size. This follows from the 
commutation relations of quantum angular momentum operators |?], [l^, [13]. 

The orientation, #(xj, t), of a Skyrmion centred at a point is described by 
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the following effective action^, |l2| |: 

^IiOfaQtfefat) - ^cos(0(x 4 ,t)-0( Xj ,t)) . (1) 

i <i,j> 

li is the moment of inertia of the i th Skyrmion and is the stiffness to relative 
rotations of neighbouring Skyrmions. The first term in Eq.(|l|) arises due to the 
change in energy of a Skyrmion when its size fluctuates; AE = I~ 1 5s 2 /8. Since 
the z-component of spin and orientation are conjugate coordinates, a cross- 
term i5s(iii, t)dtO(xi, t)/2 appears in their joint effective action. Integrating out 
5s(xi, t) gives the first term in Eq.((l|). Clearly, 1/4/ is the second derivative of 
the Skyrmion energy with respect to its spin. / is related to the Skyrmion size 
(in fact / = 24^b(7_B/s[Q) and, in the absence of disorder, is the same for all 
Skyrmions. Correlation functions involving Ss may be calculated using Eq.(|l|) by 
making the replacement £s(xj, t) — > 2ildt0(x, t)- the result of a simple Gaussian 
integration over Ss. The second term in Eq.(Q) is an effective dipole- interaction 
of Skyrmions due to the energetics of overlapping Skyrmion tails |J, [ll|. In a 
square Skyrmion lattice Jl^], is independent of lattice site. A continuum 
limit may be taken where 9i is replaced by a staggered field, 9i — ► Oi +rjiir, with 
T)i = 0, 1 on adjacent sites, and 9i — 9j is replaced by a derivative; 

S = \j ^^-)(^ 2 + ^ 2 )0(-k,^). (2) 

The frequency integral in this expression is shorthand for a Matsubara sum- 
mation at finite temperature and the momentum integral is over the Brillouin 
zone|l5|. A factor of the Skyrmion density Svp has been introduced, where 5u is 
the deviation from filling fraction v = 1 and p is the electron density. There are 
a few caveats to the use of Eq.(||). We defer discusion of these until later. This 
model is perhaps most familiar as an effective theory of the Josephson junction 
array. In this case, 9 is the phase of the superconducting order parameter and 
its conjugate coordinate is the charge of the superconducting junction. 

In order to calculate properties of the Skyrmion lattice, we must relate fluc- 
tuations in the Skyrmion orientation to fluctuations in the orientation of local 
spin. We use a coherent-state representation of the polarization of the local 
spin, via an 0(3)-vector field n(x, t). The static spin distribution at a point x 
relative to the centre of a single Skyrmion is denoted by n(x) and its in-plane 
components by n x + in y = n r e l ^° . The in-plane components of local spin at a 
point x, in response to rotational fluctuations of a Skyrmion centred at x», are 
given by 

n r e t9 {x,t) = 

We have used the conjugate relationship between Skyrmion spin and orientation 
in writing down this expression. In a distribution of many Skyrmions, one must 




n r (x 



2il 



dn r (x - 



e ie+i ^( Xi ,t). (3) 



3 



in principle sum the contributions of all Skyrmions to the fluctuation in spin at 
the point x. However, in the dilute limit in which we are performing our explicit 
calculation, Skyrmions are exponentially localized. The dominant spin fluctua- 
tions occur near to the centre of Skyrmions and so, to logarithmic accuracy, the 
local fluctuations at a point x are due only to the nearest Skyrmion. 

We will calculate the nuclear relaxation rate due to low energy quantum 
fluctuations of the Skyrmion lattice; 

JL = T 7 hm [ MlMS+MS-i-*,-*)), (4) 
Ti w-vo./ (2tt) 2 w K ' 

where 7 is the hyperfine coupling constant. Other physical observables, such 
as the temperature dependence of magnetization, (M) = J dtd 2 x(n z (x., t)), may 
be calculated similarly. Our first task is to replace the expectation of the spin 
raising and lowering operators in Eq.(Q) by correlators of the Skyrmion orien- 
tation. Substituting from Eq. (||) into Eq. (Q) and ignoring terms higher order in 
frequency, the nuclear relaxation rate at a point x is 

1 f „\_rr.. i:„ 2,„\ f d2k Zm(e td (k,u)e- l9 (~k,-iu)} 



= T7 lira in r (x) / — — - - . (5) 

This takes the form of a correlation function of the Skyrmion orientation, multi- 
plied by a profile function characteristic of the Skyrmion groundstate. The av- 
erage rate is given by integrating this over an area containing a single Skyrmion 
and multiplying by the Skyrmion density 5vp. The result is identical to Eq.(||) 
with the replacement nji(x) — > 5vn 2 . n 2 — p J d 2 xn 2 (x.) is a number character- 
istic of a single Skyrmion. For a pure Skyrmion spin distribution |lC| , n 2 — 2s to 
logarithmic accuracy in the Skyrmion spin, s = p J d 2 x(l — n 2 (x)). Notice that 
radial fluctuations of local spin contribute only to higher order in frequency and 
have been neglected in writing down Eq.(Q). 

The problem of finding the nuclear relaxation rate has now been reduced to 
evaluating the correlation function in Eq.(|J) using the effective action Eq.(^). 
This is rather tricky. The 0(2)-quantum rotor, described by Eq.(^), is quantum- 
critical, ft is necessary to employ a non-perturbative scheme, such as f/N or 
epsilon expansions, to calculate in the quantum-critical regime of this model. 
Here we use the result of the 1 /N expansion of Chubukov et q?|i^| . 

An important feature of the effective action, Eq.(j|), is that it displays a 
zero-temperature phase transition. For I J < f, the Skyrmion moment of 
inertia is sufficiently small that quantum fluctuations destroy long range or- 
der even at zero temperatures. For I J > 1, the T — groundstate has an 
infinite correlation length and is ordered. Notice that arbitrarily small tem- 
peratures destroy this long-range order even when I J > f. This arises due 
to the interplay between fluctuations and constraint and can be seen in a 
simple mean-field calculation |f6|]. In the 0(2) representation of Eq.(|^), S = 
J dtd 2 x [n(Id 2 + Jdl)n + A(x, t)(n 2 - I)] , where n(x, t) is an 0(2)-vector field 
and A(x, t) is an auxiliary field that imposes the constraint n 2 = 1. Impos- 
ing the constraint at mean-field level, (n 2 ) = 1, determines a temperature 
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Figure 1: T-dependence of the Gap. 
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Figure 2: Phase diagram of the 0(2)-vector model. 



dependent gap, A(T). The spin correlations decay exponentially on a length 
scale £(T) = \J J/X(T). The results of such a calculation are sketched in 
Fig.l. Above a temperature of about Tqc = \2nJ — i?max| (The cut-off, 
-E-max = 2iTy/~J/I corresponds to fluctuations with momentum at the Bril- 
louin zone boundary, fc m ax = ^VpdV|i5|]), the gap/correlation length devel- 
ops a universal temperature dependence. In this region, thermal and quan- 
tum fluctuations are of similar importance and are very difficult to disentan- 
gle. This crossover to universal high-temperature behaviour from distinct low- 
temperature behaviours is a feature of all correlation functions of Eq. (^) . It is 
usual to summarise this behaviour by the phase diagram sketched in Fig. 2. In 
this figure, J(T) is the renormalized stiffness in the ordered phase and A(T) is 
the gap in the paramagnetic phase. 

The correlation function required in Eq.(||) has been calculated in Ref.JT6[ 
by means of a 1/N expansion, with the result 
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Figure 3: Temperature of the nuclear relaxation rate, 1/T\. 



where 77 is a number close to zero. Substituting these results into Eq.(|5|) we 
obtain 



at high temperature. The full behaviour is sketched in Fig. 3. The kink at Tkt is 
due to the discontinuous change in spin stiffness seen at the Kosterlitz-Thouless 
transition. This effect is not seen in the 1 /N expansion of Ref . [|l6] and must be 
calculated by some other means. 

The nuclear relaxation rate obtained here is very different from that obtained 
in Ref. ]?j . As pointed out in the introduction, this is due to the unphysical limit 
T/lo — > used in Ref.0. Nevertheless, it is instructive to see how the results 
of Ref. [Q relate to the present formalism. Since T <C lu in their work, Cote et 
al consider fluctuations on length scales much less than the correlation length. 
The groundstate is ordered and a spinwave expansion may be used. Long wave- 
length fluctuations lead to a superposition of orientations- rotational averaging- 
the immediate consequence of which is that (e ie (k, uu)e~ 10 {— k, — uj)) = 0. Re- 
turning to the substitution of Eq.(||) into Eq.(^), we must retain terms to next 
order in frequency. This is a cross-term between radial and transverse fluc- 
tuations and involves a correlator (dt0e ie (k, uj)e~ ie (~ k, — u>)). Evaluating this 
correlation function via a T = spinwave expansion of the effective action, 
Eq.(||), reproduces the result of Ref. (up to a numerical factor due to our 
estimate of the Skyrmion profile function, n r dn r /ds). The zero-temperature 
phonon contribution may be calculated similarly. 

Up to now, we have assumed that the Skyrmion lattice is pinned by disorder 
and that phonons may be ignored as a consequence. The situation is rather sub- 
tle and a fuller discussion is appropriate at this juncture. A quadratic effective 
action for phonons of the Skyrmion lattice is knownj3|, 0, p|. ft is identical to 
that of an electronic Wigner crystal in a magnetic field with a vanishing effective 
mass[jl^]. Under this effective action, Skyrmions move in small ellipses with a 
frequency Wk ~ |k| 3 / 2 /S and major axes orientated transverse to the phonon 
momentum. 

At finite temperature, the occupation of transverse phonons is infra-red di- 
vergent. This divergence is restricted by interactions between phonons arising 
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from the non-harmonicity of the Skyrmion interaction. Unlike fluctuations in 
orientation, where the spinwave interaction is due to a topological constraint, 
these phonon interactions are non-universal (at best the universality is hidden 
in the details of the groundstate spin distribution and effective Skyrmion in- 
teraction potential). The resulting physics is very similar to that discussed for 
the rotation mode above; the phonon system is quantum-critical and has low- 
temperature ordered and quantum-melted^ phases and a high-temperature 
quantum-critical regime. 

Even this is not the full story. Although the low energy dispersions of 
phonons and Skyrmion rotations are independent j?]] , non- linear interactions ex- 
ist between these modes. The lattice stiffness is, in part, due to the dipole 
interaction of Skyrmions and is affected by fluctuations in orientation. Simi- 
larly, the dipole interaction between Skyrmions is strongly dependent upon the 
separation of Skyrmions and is affected by phonons. These non-linearities occur 
on the same footing as the phonon-phonon interactions and rotation-rotation 
interactions. The full quantum-critical structure of the multi-Skyrmion system 
is tremendously complicated. 

The position taken here in neglecting this wealth of structure is that the 
Skyrmion lattice is pinned by disorder and the phonon spectrum gapped. Phononic 
fluctuations are suppressed at low temperatures and the associated critical struc- 
ture occurs at higher temperature. The residual effect of phonons is a slight ther- 
mal renormalization of the rotational stiffness. The slight distortion in static 
positions of Skyrmions, in response to the disorder potential, gives a small ran- 
dom contribution to the stiffness J. This randomness produces a small region 
of Bose-glass phase at low temperatures, intervening between the paramagnet 
and renormalized classical regimes. For weak disorder, this phase only affects 
the physics very close to the critical point and does not affect our conclusions. 

We now turn to a discussion of the experimental implication of the above 
calculations. Detailed measurements of \jT\ have been carried out by Bayot et 
a?|Q. Above 40mK, Ti is independent of temperature. This is consistent with 
the rotational degrees of freedom being in their quantum-critical regime. Values 
of / and J for this system extracted from TDHFA calculations put I J very 
close to 1. The system is close to criticality and the crossover to the quantum- 
critical regime occurs at correspondingly low temperature. At AOmK there is 
an abrupt step in 1/Ti(and attendant peak in heat capacity). This is consistent 
with a Kosterlitz-Thouless transition in the orientational order (notice that the 
crossover temperature \2ttJ — £"max| may be much less than Tkt = 27rJ and 
so the behaviour may be quantum-critical either side of the transisiton) . There 
are a number of other candidate transitions^, [|, 0, however, and it is not easy 
to discriminate between them. Considerations along the lines of those presented 
here allow some elaboration, but this is necessarily rather speculative and we 
refrain from its discussion at present. We may make some firm predictions 
for nuclear relaxation measurements below AQmK . Changing the deviation in 
filling fraction or using tilted filed measurements both change the parameter I J 
and allow exploration of the phase diagram shown in Fig.l. The divergence or 
otherwise of the nuclear relaxation rate as temperature is reduced to zero should 



7 



give a clear indication of the quantum-critical structure. 
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